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Abstract 

We investigate the role that noise plays in the hysteretic dynamics of a suspended nanotube 
or a graphene sheet subject to an oscillating force. We find that not only the size but also the 
position of the hysteresis region in these systems can be controlled by noise. We also find that 
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nano- resonators act as noise rectifiers: by increasing the noise in the setup, the resonance width of 
the characteristic peak in these systems is reduced and, as a result, the quality factor is increased. 
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I. INTRODUCTION 



Recent theoretical^ and experimental woriP^ have shown that reducing the dimension 
and the size of a material with a vibrational degree of freedom enhances the role of the non- 
linearities in the dynamics. Theory and experiments have demonstrated noise squeezing^ 
and stochastic resonanc e 1 4 1 5 1 9 1 in nanomechanical resonators. These systems are used as preci- 
sion measurement devices for applications such as mass and force sensing,^ quantum motion 
detection and radiofrequency signal processing.^ Resonators consisting of suspended nan- 
otubes or graphene sheets can be used as extremely accurate mass sensors with yoctogram 
resolution.^ Models based on driven Duffing oscillators with nonlinear dissipation are used 
to interpret measurements in these devices.™ Nonlinear terms in the model are responsible 
for bistability in the oscillatory response as a function of the driving frequency and of quality 
factors dependent on the driving force. 

Motivated by experiments on nanotube and graphene resonators] 7 * 8 * 12 ! we study in this 
work the dynamics of graphene and nano tube resonators based on a suspended wired subject 
to an external oscillating driving force. We take into account the nonlinear damping term, 
the Duffing nonlinearity and we add white noise forcing not considered in the model used to 
interpret the experimental results. Given the use of these resonators resonators as accurate 
mass sensing devices, it is important to include in the dynamics all major effects in order 
to describe more accurately both, the system characteristics and the different quantities 
inferred from it, such as, mass, quality factor, dissipation coefficients, etc. 

Noise is always present in a physical system. For linear systems, noise is commonly 
regarded as having a destructive but relatively innocuous effect, blurring our view of a system 
but having no effect on the underlying processes involved. In nonlinear systems, a driving 
white noise term can modify drastically the deterministic dynamics. It can shift bifurcation 
points or induce behaviors with no deterministic counterpart, thereby affecting the dynamics 
and the quantities measured or inferred indirectly from the experimental measurements. 
Previous theoretical work showed that the main sources of noise in linear simple harmonic 
oscillators are thermomechanical noise, temperature fluctuations and adsorption-desorption 
noise.--* 1 The analysis of the sources of noise in nonlinear resonators is beyond the scope of 
this work. However we will consider two different cases. On the one hand, we shall assume 
that the sources are similar to those of linear resonators and discuss the general case of 
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additive white noise, without worrying about the causes of it. On the other hand, a local 
nonlinear damping (dissipation) means that we may need adding an appropriate white noise 
force to the equations according to the fluctuation-dissipation theorem. This will result in a 
special type of multiplicative noise. In both cases, considering the effect of noise will allow 
to estimate the values of the cubic nonlinearity with higher precision than in the noiseless 
case. Our results may also shed some light for other possible applications. 

The rest of the paper is as follows. In Section [TT] we give a general discussion of the 
characteristics of nonlinear resonators, focusing on the role that the nonlinear dissipation 



coefficient plays on the resonator dynamics. In Section III , we add an external white noise of 
fixed strength to the resonator system of Section ITT] and discuss the subsequent modifications 



of its main characteristics. In Section IV we discuss the effect produced by an internal noise 
satisfying the fluctuation-dissipation theorem. Finally, in Section [V] we present the main 
results of this work and the implications that they have in experimental measurements. 



II. THEORY 



We first discuss the properties of mechanical resonators based on a suspended nanowire 
or graphene strip that is doubly clamped and subjected to an oscillating external force 
with frequency /. Since we only care about the dynamics of the fundamental mode, we 
can consider the graphene membrane as a one dimensional object. In this regime the only 
difference between the nanowire and the graphene strip is given by the parameter values of 
each oneP 

As has already been shown experimentally,^ the dynamics of these resonators is highly 
nonlinear, as they present a force-dependent quality factor and, in some cases, hysteresis in 
the oscillation amplitude as a function of the driving frequency. The resonator dynamics 
can be described by the following equation of motion^ 

mx = —kx — ax 3 — 72 — f]x 2 x + Fdrive cos(27r ft), (1) 

where m is the effective mass, a is the Duffing term, 72; and fjx 2 x are the linear and nonlinear 
damping terms, and (k + ax 2 ) is the amplitude dependent spring's stiffness with a spring 
constant k = mu> 2 (u>q ~ 2irf is the natural frequency of the resonator fundamental mode, 
which is close to resonance). Depending on the sign of a, the Duffing term modifies the 
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stiffness of the resonator (more stiff for a > 0, softer for a < 0). Equivalently, we can work 
with the nondimensional equation 

x = — x — Q x — x 3 — i]x 2 x + Fn cos(fit) 

with ^2) 

~ / 2 + + ~ n-i T t? Fdrive n 2?r ^ 
x ^i/amwg, t = u t, Q = , r] = , F D = — 5— */— it = , 

v muj a Uq V 171 w o 

which will be used all throughout the paper. This means that all the quantities appearing 
in this work are dimensionless except when specified otherwise. 

Duffing resonators present a characteristic resonant line shape. They have a resonant peak 
for the maximum oscillation amplitude \x\ at the driving angular frequency Q = Q res and a 
force-dependent resonance width AQ. The latter is the difference between the frequencies 
of the minima flanking the resonant peak, and it is constant for linear resonators. When 
the changes to the resonance due to dissipation or noise are smaller than the resonance 
width, the standard definition of quality factor is Q = uqE/ (dE/dt) = 1.09fi/Af2 (E the 
mechanical energy at a given time and (. . .) denotes time- averaging over a time scale long 
compared with the oscillation period but sufficiently short that the decay of the amplitude 
is negligible). See the Supplementary material in Ref. EJ For linear resonators the resonant 
frequency Q res and the quality factor, Q, are independent of the driving force^, whereas 
for nonlinear resonators the energy dissipation and hence the quality factor depend on the 
oscillation amplitude. Numerical simulations show that the resonant frequency Q res increases 
with increasing r] (Figure [l^i) and Fjj (Figure [T}d). The width of the characteristic peak also 
depends on r\ and F D . 

The dependence of Q on the driving force and on the nonlinear dissipation term is not 
the only characteristic of Duffing resonators. One important feature of these systems is the 
presence of hysteresis in the oscillation amplitude as a function of the driving frequency or 
the driving forceP Figure [T] shows that the system has more than one possible oscillation 
amplitude for a range of values of Q. In the bistable region there are two stable fixed 
points corresponding to periodic orbits in the case of periodically driven wires, one of larger 
oscillation amplitude (higher in energy) and one of smaller oscillation amplitude (lower 
energy). There is also one unstable periodic orbit located between the stable ones. The 
dark line (blue line in the color version) on Figure [l] corresponds to unstable periodic orbits 
while the light (green) line corresponds to the stable ones. 




FIG. 1: (Color online) Maximum oscillation amplitude, | function of in the case of zero 

noise: (a) for a fixed value of the driving force Fd = 0.2 and different values of ij, from bottom to 
top rj = 1.2, T) = 0.5, 7] = 0.2; (b) for a fixed value of the nonlinear damping coefficient r\ = 0.5 
and different values of Fjj from bottom to top Fd = 0.05, Fd = 0.1, Fd = 0.2. Empty circles 
correspond to stable orbits and solid triangles to unstable orbits. 

The two key elements in the hysteretic behavior, discussed so far in both theory and 
experiments, are the Duffing term, responsible for the system bistability, and the nonlinear 
dissipation term fjx 2 x, which controls the size of the hysteresis region. In the limit of weak 
linear damping, perturbation theory based on the small parameter e = Q^ 1 <C 1 shows 
that it is necessary to have 77 < v^3 for the solutions of equation ^ to exhibit bistability.^ 
However this condition only provides an upper bound. Numerical simulations show that 
the maximum rj for which bistability is present, t]m, depends on the driving force amplitude 
Fd- Figures |2^i and[2]D depict tjm as a function of Fd and the size of the hysteresis window 
(the range of frequencies at which there is hysteresis), AH, as a function of r]. Note that 
77m(0.2) = 1.325, indicated i n Fig. [2^,, is appreciably smaller than the upper bound y/3, 
whereas for a much smaller Fd = 3.8 x 10~ 3 , t/m = 1.642 as indicated in Fig. [2b is closer to 
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FIG. 2: (Color online) Bistability zone as function of nonlinear dissipation amplitude rj and Q = 2irf 
for fixed values of the driving amplitude: (a) Fu = 0.2 and (Jo)Fd = 3.8 x 10~ 3 . The projection of 
the red curve on the O axis represents the minimum and maximum values of at which the system 
is bistable. The size of the hysteresis window, AH = Vt max — f2 m j n , depends on both the driving 
force and the nonlinear dissipation parameter [compare panels (a) and (b)]. For example, for 
T) = 0.5, the minimum (maximum) frequency at which hysteresis can be observed is ^l m in = 1-262 
(resp. Qmax = 1.357) in panel (a) and 0, m in = 1-02 (resp. O max = 1.034) in panel (b). As rj 
increases, AH gets reduced until the inflection point, £l ma x = ^min, at which hysteresis completely 
disappears. 

y/3. Increasing 77 reduces the frequency range at which the system is bistable. For a given 
driving force the minimum frequency Q m i n at which there is bistability does not vary 
too much, whereas the maximum frequency Q m ax at which hysteresis can be observed moves 
toward smaller values. In Figure [T^,, fl m i n = 1.257 and Q max is seemingly unbounded for 
77 = 0. For 77 = 0.5, fl m in has moved slightly to £l min = 1.251 while the maximum frequency 
has been considerably reduced to fl ma x = 1.357. This reduction continues until 77 reaches 
a critical point where Q m i n = £l m ax and bistability completely disappears. This happens at 
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r] = 1.325 in the case of F D = 0.2 (Figure^), and at rj = 1.642 for F D = 3.8 x 1(T 3 (Figure 
^p). Eichler et al estimate fj — 7.9 x 10 5 kg/(m 2 s) > a^/3/u (corresponding to i] > y/3) for 
carbon nanotubes and fj = 1.5 x 10 5 kg/(m 2 s) < ay/Z/ojQ (rj < y/3) for graphene resonatorsP 
So far the only parameter controlling the size of the hysteresis region or the dependence 
of Q on the driving force has been the nonlinear dissipation coefficient rj. This coefficient 
is an intrinsic property of each material and is fixed for every set up. However there is a 
parameter that may give rise to similar effects and that can be tuned or controlled externally. 
In the next section we explore the effect that the addition of noise has on the dynamics of 
the system and show that it has similar properties as those found for 77. 



III. EXTERNAL NOISE: RESULTS AND DISCUSSION 

According to the fluctuation-dissipation theorem, noise accompanies the process of energy 
dissipation. Thus we can expect similar behaviors as those related to 77 when noise is taken 
into account. For such purpose we add an external white noise to equation ^ 

x = —x — Q x — x 3 — i]x 2 x + Fd cos(fit) + <x£(t), (3) 

(£(*)> = o, mW)) = Kt-t'), 

which is presumably larger than the intrinsic noise. Numerical simulation of this stochastic 
equation shows that the noise shrinks the hysteresis window (Figures |3j |4^,b) and reduces 
the resonance width (Figure [5]). 

For fixed values of Fd and 77, the hysteresis window is reduced by increasing the noise 
strength, a. Figure [3] shows that, at rj — 0.5, the hysteresis window has been reduced from 
AH = 0.098 at a = to AH = 0.04 at a = 0.014. In contrast with the effect of increasing 77, 
increasing a increases fast the minimum frequency Vt min at which bistability starts, whereas 
the maximum frequency VL max is smaller than that of the deterministic case but it does not 
change than much with a; see Figure |4| Let us now compare the effects of rj and o on 
the hysteresis window. In Figure |2j we decrease the window size by increasing rj from zero 
to 7] = 0.7671 at zero noise, whereas in Figure [4] we do so by increasing a from zero to 
a = 0.014 at 77 = 0.5. In both cases AH = 0.04 is the same, however Q max and Vt min are 
different: (i) fl m ax = 1-299, Vt min = 1.259 for 77 = 0.7671 and a = 0, and (ii) fl max = 1-318, 
and flmin = 1-278 for a = 0.014 and 77 = 0.5. The hysteresis window has shifted towards 
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FIG. 3: (Color online) Curve of the maximum oscillation amplitude as a function of fi at zero 
noise, (solid line) and at a = 0.014 (doted line) for a given tj = 0.5, Fjj = 0.2, Q" 1 = 5 x 10~ 4 . 
The area under the big solid square corresponds to the bistability zone at zero noise, and the area 
under the smaller doted square represents the bistability zone for a finite noise amplitude a 



higher values of f2. The critical value at which hysteresis disappears, f2 max = $l m %ni occurs 
at = 1.251 for a = and r] = 1.325 (Figure [2]). However, for r] = 0.5 and a = 0.014, 
it occurs at Q — 1.369 far above the original zone of hysteresis at zero noise (shaded area 
on Figure [4]). Summarizing, increasing a not only reduces the hysteresis window (as it has 
already been found experimentally),^ it also shifts the window towards higher frequencies. 
This is so to the point that, just before the critical frequency fl max = ^min, hysteresis is 
found on a frequency region in which there is no bistability at zero noise (shaded area on 
Figure s|4^ and]!^). 
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FIG. 4: (Color online) Bistability zone in the plane (a, 0) for fixed nonlinear damping r] = 0.5 
and (a) Fjj = 0.2, (b) Fd = -02. The projection of the curve on the O axis gives the minimum 
and maximum frequencies at which there is hysteresis. The shaded zone corresponds to values of 
O at which there is no hysteresis for a = 0. This region gets smaller as the noise is increased. 
Comparing (a) and (b), we can observe that this region also gets smaller for smaller values of Fd 

The fact that the hysteresis window moves towards higher frequency values when a 
increases is more than just a curiosity. It is telling us a lot about what is going on in the 
system: it means that noise is feeding energy into the system. This can be understood with 
the following argument. In the hysteretic region there are three periodic orbits, one unstable 
orbit separates stable orbits with larger and smaller amplitude (therefore with higher and 
lower energy), as seen in Figures [l] and [3] When the hysteresis zone is reduced, only one 
of the stable orbits survives. The other stable orbit coalesces with the unstable one and 
both then disappear. As rj increases the surviving orbit has the smaller amplitude: as rj 
increases, the system is dissipating more energy and its dynamics selects lower energy orbits. 
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FIG. 5: (Color online) Maximum oscillation amplitude X max as a function of Q, for 7] = 1.6 and 
Fq = 0.0038. bottom solid line correspond to zero noise and top solid line a = Fo, in this regime 
hysteresis has completely disappeared. The last minimum before the curve increases monotonically 
to its global maximum are located at f2 = 1.1 for a = 0(bottom solid line), at fi = 1.055 for 
a = 0.05Fd (solid triangles), and for a = 0.5Fd (stars) and at 17 = 0.97 for a = Fd (top solid 
line). The first minimum before the curve decreases are located at £1 = 0.885 for a = (bottom 
solid line), at Q = 1.03 for a = 0.05-Fo (solid triangles), and for a = 0.5Fd (stars) and at = 1.026 
for a = Fd (top solid line) 

In contrast to this, the larger amplitude orbit survives when a increases because the noise 
is feeding energy to the syste m making more likely for the system to oscillate around the 
higher energy periodic orbit. It also turns out that the quality factor increases with noise 
which enforces this argument. 

Even more important than the reduction of AH is that increasing the noise strength 
makes the resonance width slightly thinner. In the absence of noise the two minima flanking 
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the peak are hardly perceptible. The two deep minima observed in the experiments of Eichler 
et al in Ref. [H] make us suspect that noise was present in those experiments. This effect 
that resonance width shrinks with noise is always present, even in systems not exhibiting 
hysteresis. When a increases, the two minima flanking the peak become closer and deeper. 
On Figure [5j we see how the two minima are almost imperceptible and are located at 
Q = 0.88 and Q = 1.05375 at zero noise (green crosses). For a = Fd, the two minima are 
more pronounced and are located at — 0.9765 and Q = 1.03125. The resonance width 
has been reduced from AQ = 0.173 at a = to AQ = 0.05 at a = Fp (blue stars). This 
shrinking can be se en in Figure [5] for more values of a. This indicates that the presence of 
noise is affecting the value of the quality factor Q. Since the reduction/enlargement due to 
noise and dissipation are smaller than the resonance width, the quality factor is^ 

« = n »^" 1 - M M- (4) 

where E is the mechanical energy at a given time t. According to the Supplementary 
material in Ref. El the approximation Q « 1.09fi/Afi holds in the limit of very large quality 
factor. The noise reduces AQ and therefore it enhances the quality factor. This contrasts 
with the behavior of linear systems. In them, external noise increases energy dissipation and 
makes these systems less efficient. However for our nonlinear system, the random energy 
provided by noise gets absorbed by the system and converted in useful mechanical energy, 
thereby yielding a higher quality factor. In this sense, the nonlinearly damped driven Duffing 
oscillator is a noise rectifier that improves its performance by using energy provided by the 
external noise. 

The reduction of the frequency width is also telling us that the value of r] calculated from 
the experiments of Eichler et ai s (in which noise is certainly present) may be smaller than 
the real one, i.e., r]\ u=0 > r]\ a>0 . The value of r\ in Ref. [8] is calculated from the expression 

mA/\ 3 (F D S 2 



where the dependence of AQ on a has not been taken into account. We can include this 
dependence on equation ^ and let i](a) be given by ^ for a fixed value of a > 0. We can 
rewrite AQ = AQ \ a=0 //(cr), with f(a) = (AQ \ a=0 /AD, \ a> o). The factor f(a) indicates 
how the resonance width changes with noise: /(0) = 1 and f(a) > 1 for a > 0. Using ([5]), 
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FIG. 6: (Color online) Maximum oscillation amplitude X max as a function of 0, for r/ = 1.6 
and Fd = 0.0038 when the noise strength is given by the fluctuation-dissipation theorem, a 2 = 
25(r/x 2 + 1/Q), with 5 = ak B T/(k 2 ). From bottom to top T = 0, T = 50K, T = 100K, T = 200K, 
T = 300K. It can be observed that the resonance width is reduced as T (hence also a) increases, 
although in this case it is more difficult to identify the two minima flanking the peak. 

we see that 

rtip) = J7%< V(0). (6) 

This means that the estimated value i](a) is smaller than the real nonlinear dissipation 
coefficient of the deterministic equation, rj = 1](0). The real value of rj is screened by the 
effect of noise. 
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IV. INTERNAL NOISE 



So far we have discussed a general case of external white noise, using a constant noise 
amplitude a, without specifying the nature of it. However, in our system there are different 
sources of dissipation that contribute to 7 and fj in 0. Let us assume that the dominant 
noise is thermo- mechanical.^ According to the fluctuation-dissipation theorem, there should 
be a stochastic white noise forcing term in ^ with an appropriate strength, cr 2 = 25(r]x 2 + 
l/Q), where 5 = akBT/k 2 . We have performed numerical simulations of ([3j with this 
internal noise for a given value of 77. The results obtained are qualitatively similar to those 
discussed above. Figure [6] shows that increasing the temperature reduces the resonance 
width (similar to the effect of increasing the fixed strength a of an external noise), thereby 
increasing the quality factor. Similarly, the hyste resis frequency range diminishes with all 
the consequences already mentioned. 

The measurements of the quality factor made by Eichler et al are not affected by these 
results, since they calculate Q = 1.090o/Afi and take AQ from the experimental mea- 
surements which are affected by the existing noise. However the inferred value of fj may be 
affected by including noise. The measurements in Ref. [8]were interpreted assuming that only 
the nonlinear dissipation fj was responsible for the shrinking or expansion of the resonant 
peak. Since noise also has an effect, the value of fj proposed in Ref. could be smaller than 
the real one. 

V. CONCLUSIONS 

We have found that the presence of noise modifies in a qualitative way the dynamics of 
the nonlinear Duffing resonators, and it gives rise to some effects similar to those attributed 
to the nonlinear dissipation coefficient. Noise not only reduces the window of hysteresis but 
also shifts it to higher values of the driving frequency. Since the value of rj is inherent to 
the material and can not be tuned, there is not much room for modifying the quality factor 
rather than increasing Fp which cannot always be done experimentally. However controlling 
the sources of noise allows to increase or decrease the bistability regions and even to shift 
them to higher frequencies within easy experimental reach. 

We have also found that the resonant width can be controlled by noise and hence the 
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quality factor can be increased by increasing the sources of noise in the system. This means 
not only that higher quality factors can be reached by increasing the sources of noise, but 
also, and even more important, it means that nanoresonators act as noise rectifiers. This 
is the main result of this work, to show that nonlinear nanoresonators are able to convert 
random energy from the environment into useful mechanical energy. 

Even though there are different sources of noise present in these systemsj 13 * 15 4n this work 
we have discussed mostly the case of white noise with constant strength a, leaving its nature 
unspecified. However, we also performed numerical simulations assuming that the dominant 
noise is of thermo-mechanical origin and satisfies the fluctuation-dissipation theorem (the 
noise strength a = \j25{rjx 2 + 1/Q), with 5 = aksT/k 2 , depends on the temperature and 
on the oscillation amplitude). The results are qualitatively similar to those for constant 
strength white noise. 
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